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I .   INTRODUCTION 

A .   BACKGROUND 

Fleet  Numerical  Oceanography  Center  (FNOC)  provides 
numerical  and  oceanographic  products  for  use  by  subordinate 
and  individual  commands.  FNOC  is  the  master  computer  center 
for  the  Naval  Oceanography  and  Meteorological  Support  System 
(NOMSS) .  Operationally,  FNOC  falls  under  the  Naval 
Oceanography  Command.  Data  is  received  from  around  the  world 
and  is  used  to  produce  a  wide  variety  of  products  designed  to 
benefit  the  military,  particularly  the  U.  S.  Navy.  Several 
broad  classifications  of  the  types  of  products  exist.  These 
are  atmospheric  weather  conditions,  ocean  weather  conditions, 
radar  propagation  data,  and  underwater  conditions. 

Atmospheric  weather  condition  data  is  used  primarily  by 
the  Optimum  Path  Aircraft  Routing  System  (OPARS)  .  This  system 
is  a  computer  model  that  receives  requests  for  flight  plans 
directly  from  individual  users.  The  program  then  processes 
the  request  and  sends  out  an  optimal  flight  plan  based  on 
actual  or  climatological  weather  conditions.  The  choice  of 
actual  or  climatological  conditions  is  based  on  the  requested 
date  for  the  flight  plan.  FNOC  estimates  that  this  computer 
model  saves  the  military  ten  million  dollars  annually  in  fuel 
costs . 


Ocean  weather  data  is  provided  to  two  centers,  one  in 
Norfolk,  Virginia  and  the  other  in  Pearl  Harbor,  Hawaii. 
These  centers  use  this  data  to  provide  USN,  USNS,  and 
contracted  vessels  with  Optimum  Track  Ship  Routes  (OTSR) . 
Currently,  these  routes  are  manually  generated  and 
distributed.  In  this  case  FNOC  estimates  that  the  annual 
savings  in  fuel  and  damage  costs  is  seventeen  million  dollars. 

Radar  propagation  data  is  computed  using  atmospheric 
weather  conditions.  The  daily  radar  range  estimates  are 
generated  by  onsite  computer  models  and  distributed  to 
reguesting  commands.  Because  of  the  diversity  of  radars  in 
use  by  the  military,  this  model  incorporates  parameters 
specific  to  the  requesting  command's  radar.  FNOC  has  no 
estimate  of  the  benefit  provided  by  these  models. 

Underwater  conditions  are  also  predicted  by  computer 
models.  "The  output  from  these  models  is  then  used  by  other 
computer  models  to  determine  estimates  of  sonar  ranges  for 
various  ship  configurations.  Ranges  are  predicted  for  both 
active  and  passive  sonars  and  sonobouys .  As  with  the  radar 
propagation  data,  requesting  commands  provide  the  type  of 
sonar  or  sonobouy  along  with  other  operating  parameters. 

Of  the  products  discussed  above,  two  are  most  beneficial 
during  peacetime.  The  first  is  the  Optimum  Track  Ship  Routing 
(OTSR)  and  the  second  is  the  Optimum  Path  Aircraft  Routing 
System  (OPARS)  .  As  FNOC  has  estimated,  these  two  products 
generate  a  savings  to  the  military  in  excess  of  seven  million 


dollars  over  FNOC ' s  operating  budget.  These  estimates  are 
based  on  conjecture  and  do  not  have  an  underlying  model  to 
support  them. 

B.   PROBLEM  STATEMENT 

It  is  desired  to  produce  a  model  that  can  be  used  to 
determine  the  benefit  gained  from  the  products  generated  by 
FNOC.  Since  OTSR  and  OPARS  appear  to  provide  the  large  bulk 
of  the  peacetime  savings,  they  will  be  addressed  in  the 
following  analysis. 

OTSR  provides  a  recommendation  to  ship  captains  and 
masters  on  the  track  that  would  save  the  most  time  and  fuel 
and  result  in  the  least  damage  to  the  ship  from  weather. 
Without  this  system,  ships  would  take  routes  based  on 
historical  climatological  data.  In  fact,  while  the  first  few 
days  of  the  initial  OTSR  uses  actual  weather  forecasts,  the 
remainder  of  the  voyage  is  based  solely  on  climatological 
data,  just  as  any  captain  would  plan  his  route.  It  is  the 
routing  updates  generated  by  the  centers  using  ongoing  weather 
forecasts  that  generate  the  most  cost  savings. 

Several  studies  have  looked  at  this  system,  the  latest 
being  completed  in  1976  by  Lulejian  &  Associates,  Inc. 
[Reference  1] .  Although  detailed,  this  study  looked  only  at 
the  costs  associated  with  the  weather  centers  that  actually 
produced  the  OTSR  and  not  the  costs  incurred  by  FNOC  in 
providing  the  required  information.   Realizing  that  there  is 


also  a  cost  incurred  in  gathering  the  information  and 
providing  it  to  the  ships,  aircraft,  and  personnel,  this  paper 
will  address  only  those  costs  incurred  by  FNOC  and  the  two 
centers.  The  reason  for  this  is  that  FNOC  is  only  one  of 
several  weather  data  collection  agencies;  the  National  Weather 
Service  and  the  Air  Force  Weather  Center  receive  the  same 
information.  Additionally,  much  of  the  data  used  by  FNOC  is 
collected  in  conjunction  with  routine  military  operations, 
with  the  exception  of  hurricane /typhoon  locator  flights. 
These  flights,  in  the  absence  of  FNOC,  would  be  conducted  for 
the  National  Weather  Service  to  provide  early  warning  to 
coastal  regions  that  may  be  affected  by  the  storm.  It  is 
therefore  concluded  that  data  collection  is  not  unique  to  FNOC 
and  will  not  be  considered  as  a  cost. 

In  order  to  determine  the  benefit  received  from  the  OTSR 
system,  it  will  be  necessary  to  determine  how  ships  would  be 
routed  in  the  absence  of  OTSR  and  how  they  would  be  routed 
with  perfect  information.  Routings  can  be  made  in  three 
different  ways.  The  great  circle  route  is  the  shortest 
distance  that  can  be  travelled,  and  is  also  the  easiest  and 
least  costly  to  calculate.  This  choice,  however,  could  result 
in  severe  weather  encounters,  thereby  negating  any  fuel 
savings  with  damage  costs.  This  route  would  provide  an  upper 
bound  on  the  cost  in  fuel  and  damage,  because  even  without 
OTSR,  a  better  route  could  be  chosen. 


The  second  alternative  is  a  route  based  solely  on 
climatology.  Routings  of  this  sort  have  been  conducted  for 
centuries.  Although  more  costly  than  a  great  circle  route  to 
compute,  nearly  any  captain  with  access  to  pilot  charts  can 
compute  a  climatological  course.  Again,  however,  the  danger 
of  encountering  severe  weather  still  exists  since 
climatological  routes  are  only  an  expectation  of  future 
weather  in  a  region. 

The  most  desirable  choice  is  a  route  based  on  perfect 
information.  If  exact  weather  conditions  could  be  predicted, 
the  optimum  route  in  terms  of  fuel  savings  and  damage 
avoidance  could  be  chosen.  The  OTSR  system  provides  a  route 
that  represents  the  cost  of  expected  damage  and  fuel 
consumption  that  lies  somewhere  between  the  cost  incurred  by 
climatological  routes  and  that  of  a  route  based  on  perfect 
weather  information.  It  is  not  the  aim  of  this  thesis  to 
perfectly  predict  weather  conditions,  but  to  determine  the 
savings  of  the  current  system  over  the  use  of  climatological 
routes.  Climatological  routes  in  this  case  provide  a  worst 
case  situation  in  determining  the  cost  of  expected  damage  and 
fuel  consumption.  Routes  costing  more  in  damage  could  be 
chosen,  but  this  is  unlikely.  Traditionally,  ships  have 
followed  established  climatological  routes  to  obtain  the  least 
cost  due  to  damage  in  the  long  run.  By  determining  the 
expected  cost  when  using  climatological  routes,  an  estimation 
of  the  benefit  gained  by  OTSR  can  be  determined.   The  major 


benefit  of  OTSR  over  climatological  routes  is  that  OTSR  takes 
into  account  current  and  forecasted  weather  in  order  to 
determine  a  route.  Additionally,  as  forecasting  continues 
during  the  route,  adjustments  can  be  made  to  take  advantage  of 
unexpected  fair  weather  in  a  region  that  would  otherwise  be 
avoided  by  climatological  routes.  It  is  the  dynamic  nature  of 
OTSR  that  allows  it  to  make  great  gains  in  damage  avoidance 
and  fuel  savings . 

OPARS  is  a  computer  program  that  provides  direct  access  to 
users  on  optimum  paths  for  aircraft  based  on  the  following: 

•  aircraft  performance  parameters 

•  weather  conditions 

•  minimum  fuel  consumption  or  least  time  enroute  for  the 
flight  requested. 

The  majority  of  the  flight  plans  are  generated  for  Navy  and 
Coast  Guard  units,  with  the  Air  Force  and  Army  making  up  about 
20  percent  of  the  requests. 

The  major  cost  savings  associated  with  OPARS  are  fuel, 
damage  avoidance,  flight  time,  and  flight  planners'  time.  As 
with  the  OTSR,  these  flight  plans  are  only  recommendations. 
Other  operational  considerations  may  preclude  the  use  of  the 
optimal  flight  plan. 

Flight  plans  can  be  calculated  in  the  same  manner  as  the 
ship  routes.  That  is,  by  great  circle,  climatological  route, 
or  perfect  information.   Currently  flight  plans  produced  by 


FNOC  are  better  than  climatological  routes  but  fall  short  of 
the  optimum  that  could  be  obtained  with  perfect  information. 
As  with  the  OTSR,  the  problem  is  to  develop  a  model  that  will 
simulate  the  route  that  would  be  chosen  if  the  OPARS  model 
were  not  available.  This  will  be  accomplished  by  modifying 
actual  flight  plans  to  determine  the  range  of  fuel  consumption 
around  the  optimum  flight  plan  chosen  by  OPARS. 

The  unaided  flight  planner  would  be  required  to  sift 
through  all  applicable  weather  information  to  determine  the 
optimum  route  by  hand.  This  would  be  the  same  information 
that  is  currently  provided  to  OPARS.  It  can  be  expected  that 
an  experienced  flight  planner  would  choose  a  route  that  is 
close  to  the  optimum  chosen  by  OPARS.  The  modified  flight 
plans  will  provide  the  range  of  fuel  usage  around  the  optimum. 
Therefore,  the  amount  of  fuel  and  the  amount  of  flight 
planning" time  that  is  saved  by  OPARS  is  a  significant  measure 
of  its  effectiveness  and  worth. 


II .   METHODOLOGY 

A.   ASSUMPTIONS 

This  section  is  a  brief  description  of  the  assumptions 
necessary  for  the  model  formulation.  A  more  detailed 
description  of  the  assumptions  made  here  follows  in  later 
sections . 

1.   Optimum  Track  Ship  Routing  System 

The  following  assumptions  are  necessary  regarding  the 
analysis  of  the  OTSR  system. 


Perfect  information  results  in  minimal  transit  cost  due  to 
damage  and  fuel  consumption. 

Climatological  routing  gives  an  upper  bound  on  cost  due  to 
damage . 

Ocean  grid  gives  rise  to  a  sparse  graph  since  all  points 
are  not  directly  accessible  from  a  given  point. 

Spruance  class  destroyer  as  a  representative  ship  for 
model . 

Wave  height  is  the  only  significant  parameter  involved  in 
ship  damage . 

Wave  heights  at  grid  points  are  independent. 

The  conditional  probability  of  damage  given  sea  height  is 
known . 

95   percent   of   routings   are   accepted  by  Commanding 
Officers . 


Further  explanation  of  these  assumptions  is  in  the  following 
section. 


2 .   Optimum  Path  Aircraft  Routing  System 

The  assumptions  necessary  for  OPARS  are  as  follows: 

•  Fuel  and  flight  planner's  time  provide  the  savings. 

•  OPARS  route  is  optimum. 

•  The  bulk  of  the  fuel  savings  is  realized  by  a  small  number 
of  aircraft  that  fly  the  majority  of  the  flight  plans. 

•  The  unaided  flight  planner  would  choose  a  flight  plan 
within  4000  feet  of  the  optimum  altitude,  normally 
distributed  about  the  optimum  altitude. 

•  OPARS  is  capable  of  calculating  the  .fuel  required  for 
alternate  non-optimal  flight  plans. 

•  All  flights  are  flown  using  an  OPARS  route. 

•  The  FNOC  weather  model  is  accurate. 

Further  explanation  of  these  assumptions  is  in  the  following 
section. 

B.   DISCUSSION  OF  ASSUMPTIONS  AND  MODELING 
1.   Optimum  Track  Ship  Routing  System 

The  difficulties  in  routing  ships  optimally  on  long 
voyages  are  numerous.  Forecasting  techniques  are  only 
accurate  for  up  to  a  few  days  from  the  forecast  date.  Since 
ships  move  relatively  slowly,  great  care"  must  be  taken  to 
avoid  placing  a  ship  in  a  situation  from  which  it  cannot 
easily  escape.  Routings  must  be  closely  monitored  and  updated 
continuously  as  weather  conditions  change.  In  this  way,  OTSR 
uses  a  somewhat  greedy  strategy  in  that  the  initial  three  to 
four  days  transit  is  based  on  forecasted  weather  and  the 


remainder  of  the  initial  routing  is  based  on  climatological 
data.  It  is  easy  to  see  that  this  will  not  always  lead  to  the 
optimum  routing  in  terms  of  fuel  and  damage  avoidance  as  could 
be  expected  with  perfect  weather  information.  By  the  choice 
of  the  initial  days  of  the  route,  future  options  to  take 
advantage  of  fair  weather  may  not  be  available.  In  other 
words,  some  damage  can  be  expected  even  on  an  OTSR  route. 
This  leads  to  the  first  assumption,  that  perfect  information 
would  result  in  minimal  cost  due  to  damage  and  fuel 
consumption  and  that  climatological  routing  would  result  in  an 
upper  bound  on  cost  since  this  is  the  worst  routing  that  could 
be  chosen  using  all  available  information  with  the  exception 
of  OTSR.  Prior  to  OTSR,  climatological  routing  was  the  best 
available  choice. 

A  second  assumption  that  is  necessary  for  the 
forthcoming  model  is  that  not  all  points  in  the  ocean  are 
accessible  from  the  ship's  current  position.  The  model  used 
to  calculate  the  cost  of  a  route  uses  a  grid  in  which  each 
point  is  separated  by  five  degrees  of  latitude  and  longitude. 
From  the  ship's  current  point,  three  points  are  accessible  to 
it,  namely  five  degrees  of  longitude  further  along  its  track 
and  its  current  latitude  plus  or  minus  five  degrees.  This 
gives  rise  to  a  sparse  graph  that  can  be  used  in  the  shortest 
path  algorithm  to  be  described  later. 

The  model  will  also  use  a  Spruance  class  destroyer  as 
the  platform  to  calculate  fuel  usage  for  the  route.   This  is 
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soon  to  be  the  most  prevalent  engineering  plant  in  the  fleet, 
and  its  fuel  usage  is  representative  of  the  fleet. 

The  only  parameter  to  be  used  in  calculating  the 
damage  that  a  vessel  encounters  on  the  route  will  be  wave 
height.  Previous  studies  [Reference  1  and  2]  have  shown  that 
wave  height  has  the  most  significant  effect  on  the  damage  to 
a  vessel.  Wave  period  and  direction  also  play  an  important 
role,  especially  if  the  period  is  such  to  cause  resonance  at 
the  current  speed.  This  problem  can  easily  be  solved  by  an 
adjustment  of  ship  speed.  Nagel  [Reference  2]  has  shown  that 
the  effect  of  this  speed  decrease  is  small  between  the  optimum 
route  and  the  climatological  route.  That  is  to  say,  it  is 
felt  that  the  benefit  of  optimum  routing  is  greater  in  terms 
of  damage  avoidance  than  in  terms  of  time  saved. 

The  probability  of  a  particular  wave  height  at  a  given 
point  oh  the  grid  described  earlier  is  derived  from 
climatological  charts  in  the  Defense  Mapping  Agency's  Sailing 
Directions  [Reference  3].  The  probability  of  sea  height  in 
these  charts  is  based  on  observed  wave  height  during  a 
specific  month  in  the  case  of  the  North  Atlantic  Ocean  or  a 
specific  season  in  the  case  of  the  North  Pacific  Ocean.  The 
use  of  a  five  degree  grid  was  chosen  to  gain  independence  of 
the  wave  heights  from  one  grid  point  to  the  next.  Allowing 
roughly  300  nautical  miles  between  points,  creates  a  large 
enough  fetch  for  seas  to  fully  develop  in  that  region  and  not 
necessarily  be  influenced  by  an  adjacent  region.  That  is,  the 
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area  of  ocean  covered  by  each  grid  point  is  large  enough  so  as 
to  maintain  its  own  sea  height  without  regard  to  adjacent 
conditions . 

Finally,  values  for  the  probability  of  damage  given 
sea  heights  for  particular  ship  types  is  not  known,  nor  is  it 
necessary  in  order  to  develop  a  relative  cost  for 
climatological  routes  over  routes  based  on  perfect 
information.  Aggregate  values  for  the  conditional  probability 
of  damage  given  sea  height  have  been  determined  [Reference  1]  . 
These  conditional  probabilities  are  based  on  historical  data 
from  July  1969  to  June  1975  from  records  of  the  Naval  Safety 
Center,  and  are  recreated  in  Table  I. 


Table  I       CONDITIONAL  PROBABILITY  OF  DAMAGE  GIVEN  SEA 
HEIGHT  WITH  AVERAGE  COST  OF  DAMAGE  [REFERENCE  1] 


Sea  Height,  x 
(in  feet) 

Conditional 
Probability  of 
Damage  Given 
Sea  Height 

Average  Damage 
per  Incident 
(dollars) 

0<x<4 

0.0000 

0 

4<x<8 

0.0001 

48427 

8<x<12 

0.0008 

12<x<16 

0.0009 

129969 

16<x<20 

0.0118 

20<x<24 

0.029 

312196 

24<x<2  8 

0.0700 

x>28 

0.2860 

340771 
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All  costs  have  been  converted  to  1992  dollars  using 
six  percent  inflation.  The  dollar  amounts  are  estimates  based 
on  the  Commanding  Officer's  assessment  of  the  damage  and  are 
therefore  not  actual  cost  to  repair  the  damage.  For  this 
reason,  it  is  felt  that  these  values  are  lower  than  the  actual 
cost  of  repairs.  Data  for  the  average  damage  per  incident  was 
available  only  for  eight  foot  increments  from  4  to  28  feet. 
The  value  in  each  eight  foot  increment  will  be  used  with  the 
conditional  probability  of  damage  given  sea  height  values 
within  this  increment.  The  probability  of  damage  given  here 
is  without  regard  to  ship  type  or  class,  but  is  pooled  from 
available  data  from  the  Naval  Safety  Center. 

As  mentioned  earlier,  no  follow-up  by  the  Naval  Safety 
Center  is  conducted  to  determine  actual  costs  of  damage  by 
unfavorable  weather  conditions.  These  figures  represent  an 
estimate"  of  the  cost  to  repair  damage.  They  do  not  include 
the  cost  in  loss  of  availability  of  the  ship's  services.  Loss 
of  availability  may  or  may  not  be  applicable.  Much  of  the 
damage  caused  by  adverse  weather  is  not  of  a  serious  enough 
nature  to  require  the  ship  to  be  taken  out  of  action  to 
repair.  The  bulk  of  the  damage  can  be  repaired  during 
scheduled  maintenance  periods  and  would  therefore  not  impinge 
on  ship  operations.  Table  II  is  a  summary  of  damage  sustained 
by  USN  and  USNS  vessels  for  the  period  from  January  1982  to 
May  1992. 
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Table  II      SUMMARY  OF  ACTUAL  DAMAGE  SUSTAINED  JANUARY 
1982  TO  MAY  19  92,  FROM  NAVAL  SAFETY  CENTER  DATA. 


Totals 

Ave .  per  Year 

Damage  Cost 

38.659  million 

3 .741  million 

Ship  Days  Lost 

172 

16.6 

Lost  Work  Days 

617 

59.7 

Fatalities 

14 

1.4 

No.  of  Incidents 

279 

27 

Figure  1  shows  the  monthly  mean  of  the  damage  for  the 
same  time  period.  The  means  follow  closely  what  would  be 
expected  during  the  winter  and  summer  seasons,  with  the 
exception  of  April  and  October.  This  was  due  to  a  single 
unusually  high  cost  in  each  of  these  months.  If  this  value  is 
eliminated  in  each  month,  the  means  are  as  depicted  in  Figure 
2. 

Table  II  and  Figures  1  and  2  are  introduced  as 
indicators  of  the  damage  sustained  by  vessels  even  while 
operating  under  the  OTSR  system.  Although  OTSR  will  be  shown 
to  be  very  beneficial,  we  cannot  assume  that  the  system  is 
perfect  and  not  without  limitations.  Even  under  a  routing 
system  like  OTSR,  some  damage  will  occur. 

OTSR  routings  are  advisory  in  nature.  There  currently 
are  no  requirements  for  Commanders  and  Commanding  Officers  to 
follow  these  routes.   In  the  case  of  OTSR,  routings  are  viewed 
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Mean  Damage  Cost  by  Month 
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Figure  2     Mean  Damage  Cost  by  Month,  Abnormal 
Values  Removed 
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as  extremely  beneficial  by  Commanding  Officers.  This  is 
evident  in  a  high  acceptance  rate  of  the  recommendations  .  An 
acceptance  rate  of  95  percent  has  been  shown  [Reference  1] . 
The  routings  are  not  perfect  though.  In  this  same  study,  it 
was  shown  that  approximately  11  percent  of  the  routed  ships 
received  routing  changes  during  their  voyages.  These  changes 
were  to: 

•  avoid  adverse  weather,  and 

•  take  advantage  of  unexpected  favorable  weather. 

The  results  of  these  route  changes  are  shown  in  Table  III. 


Table  III     SEA   STATE   ALONG 
DIVERTED  ROUTE 


Better 

82  percent 

Equal 

12  percent 

Worse 

6  percent 

As  shown,  in  Table  III,  upon  analysis  of  the  weather 
conditions  of  the  route  taken  and  the  recommended  route,  94 
percent  of  the  ships  that  were  rerouted  experienced  seas  of 
equal  or  lesser  severity.  Only  six  percent  encountered  more 
damaging  seas.  As  the  above  study  [Reference  1]  went  on  to 
show,  four  of  those  ships  rerouted  chose  not  to  follow  the 
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recommended  course  changes.  Of  these,  two  encountered  rough 
seas  and  one  ship  suffered  damage. 

An  acceptance  rate  of  0.95  will  be  used  throughout 
this  analysis.  Additionally,  ships  that  do  not  accept  OTSR 
recommendations  will  encounter  heavier  seas  at  the  rate  of  0.5 
from  observations  [Reference  1]  of  ships  that  chose  not  to 
accept  rerouting.  Those  ships  that  do  follow  OTSR  rerouting 
directions  will  experience  heavier  seas  at  the  lower  rate  of 
0.06,  reflecting  the  error  rate  in  OTSR  rerouting. 

The  following  discussion  will  describe  a  method  to 
determine  the  savings  by  optimally  routing  ships  vice  routing 
by  climatological  data. 

In  the  absence  of  OTSR,  ships  would  be  routed  with 
climatological  and  short  range  weather  predictions. 
Historical  wave  height  information  is  available  from 
climatological  charts  [Reference  3].  In  this  model,  each  five 
degrees  of  latitude  and  longitude  in  the  North  Atlantic  and 
North  Pacific  is  assigned  a  probability  of  a  specific  wave 
height.  In  the  North  Atlantic  this  is  the  probability  of  seas 
greater  than  12  feet  and  greater  than  2  0  feet.  For  the  North 
Pacific  it  is  the  probability  of  seas  greater  than  8  feet  and 
greater  than  12  feet.  Data  for  the  month  of  January  in  the 
North  Atlantic  will  be  used  to  estimate  the  savings  gained 
from  using  OPARS.  The  region  used  extended  from  08  0 °W  to 
005 °W  and  from  20°N  to  65 °N.  Data  points  were  numbered  from 
one  at  080 °W,  20 °N,  north  to  ten  and  then  starting  again  at 
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075°W,  20°N  with  eleven  and  continuing  in  that  fashion  to  160 
at  005°W,  65°N. 

As  described  earlier  the  graph  that  is  developed  is 
sparse  with  only  a  limited  number  of  points  accessible  from 
the  parent  node.  All  arcs  in  the  graph  are  directed,  and 
movement  is  allowed  only  in  one  direction.  Additionally,  all 
costs  of  moving  from  one  node  to  the  next  are  positive  so  that 
no  negative  cycles  can  develop.  To  determine  the  shortest 
path,  i.e.,  the  optimal  route  through  this  directed  acyclic 
graph,  Dijkstra's  algorithm  [Reference  4]  is  used.  The 
complete  algorithm  is  contained  in  Appendix  A.  At  each  node, 
a  random  number  is  drawn  from  a  uniform  distribution  to 
determine  the  wave  height  at  that  node.  The  climatological 
charts  provide  the  probability  of  wave  height  for  all  ocean 
areas.  In  this  way,  a  situation  is  produced  in  which  perfect 
weather  information  is  known  for  the  entire  region.  A  second 
random  number  will  be  drawn  and  compared  to  the  probability  of 
damage  given  sea  height  from  Table  I .  The  damage  and  the 
fuel  used  to  arrive  at  that  node  provide  the  cost  of  transit 
to  that  node.  The  distance  between  nodes  is  computed  using 
the  following  formula: 

£>ist=60arccos  [sinLdsinLd+cosLscosZ,dcos  ikd-Xs)  ] 

where : 


•  Ls  is  the  latitude  of  start 
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•  Ld  is  the  latitude  of  destination 

•  Xs   is  the  longitude  of  start 

•  Xd   is  the  longitude  of  destination,  and 

•  Dist  is  the  great  circle  distance  between  start  and 
destination.  [Reference  5] 

Dijkstra's  algorithm  adds  each  node  to  a  set  of 
completed  nodes  one  at  a  time.  As  the  node- is  added,  all  arcs 
leaving  the  node  are  examined  and  adjacent  nodes  are  updated 
if  an  improvement  in  the  distance  to  that  node  is  found.  The 
algorithm  computes  the  distance  from  the  source  to  all  other 
nodes.  By  choosing  the  destination  and  tracing  back  through 
the  parents,  the  shortest  path  is  defined. 

In  order  to  determine  the  cost  of  routing  ships  in  the 
absence  of  OTSR,  a  second  algorithm  is  used.  Utilizing  the 
same  random  number  seed  as  above,  wave  heights  and  the 
probability  of  damage  conditions  were  duplicated.  This  time 
though,  instead  of  determining  the  optimum  route,  a  greedy 
strategy  was  employed.  The  complete  algorithm  is  contained  in 
Appendix  B.  Prior  to  the  start  of  the  algorithm,  all  nodes 
that  can  lead  to  the  destination  are  marked.  From  the  start 
node  then,  the  least  cost  route  to  the  next  accessible  marked 
node  is  chosen.  From  this  node  then,  the  next  least  cost  arc 
is  chosen  until  the  destination  is  reached.  In  this  way,  a 
cost  can  be  determined  for  a  route  in  which  a  ship  has  placed 
itself  in  a  situation  where  high  seas  must  be  encountered  to 
reach  the  destination. 
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After  multiple  runs  of  each  algorithm,  the  cost  of  the 
damage  incurred  and  the  fuel  used  on  the  shortest  path  is 
compared  to  the  cost  of  the  non-optimal  route  determined  in 
the  second  algorithm.  A  single  route  will  be  considered.  To 
determine  an  annual  savings  it  will  be  necessary  to  multiply 
the  mean  savings  from  the  algorithm  by  99  routings  per  month 
times  twelve  months,  for  the  North  Atlantic,  and  12  6  routings 
per  month  times  twelve  months  times  2 . 5  to  compensate  for  the 
greater  distance  traveled  in  the  North  Pacific.  This  will 
provide  a  mean  savings  if  perfect  information  were  available. 
From  Table  II  it  is  known  that  3.741  million  dollars  in  damage 
is  sustained  each  year  under  OTSR.  A  yearly  mean  was 
determined  from  the  124  months  of  data  obtained.  The 
difference  between  this  value  and  the  cost  of  damage  incurred 
on  the  optimal  routes  will  be  subtracted  from  the  final  value 
determined  for  the  savings.  This  will  provide  an  estimate  for 
the  savings  possible  under  OTSR. 

2.   Optimum  Path  Aircraft  Routing  System 

In  creating  a  model  to  establish  the  benefit  gained 
from  using  a  computer  model  to  optimally  route  aircraft  over 
routing  each  aircraft  by  hand,  it  must  be  established  what  it 
is  that  provides  the  greatest  gain.  Overwhelmingly,  the 
answer  is  fuel  savings.  Unlike  ships,  aircraft  are  capable  of 
quickly  changing  course  to  avoid  adverse  weather  conditions. 
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Additionally,  routings  usually  are  much  shorter  in  duration 
and  are  therefore  able  to  take  full  advantage  of  short  range 
weather  forecasts. 

The  Naval  Safety  Center  does  not  classify  weather  as 
a  cause  for  aircraft  damage.  This  is  because  the  aircraft  is 
either  all-weather  or  it  does  not  fly  when  forecasts  show  that 
the  aircraft  would  encounter  adverse  weather.  It  is  concluded 
that  fuel  savings  from  the  optimal  routing  of  the  aircraft  is 
by  far  a  greater  indication  of  its  benefit  than  is  damage  from 
adverse  weather. 

Secondly,  it  is  assumed  that  the  route  provided  by 
OPARS  is  indeed  the  optimum  route,  since  it  is  not  feasible  to 
verify  this  by  hindsight  routing  of  the  aircraft.  In  any 
case,  OPARS  is  the  closest  routing  system  available  to  the 
true  optimum.  The  desire  is  to  determine  how  much  this  system 
benefits"  the  military  when  compared  to  the  next  best 
alternative,  i.e.,  manually  computed  routes,  and  not  to 
determine  how  much  could  be  saved  with  some  other  system. 

The  OPARS  database  is  capable  of  providing  routings 
for  over  ninety  aircraft  types.  A  limited  number  of  these 
aircraft  types  were  chosen  for  study.  Eleven  aircraft  and 
their  variants  were  chosen  for  use  in  the  analysis.  These 
eleven  aircraft  cover  twenty-one  of  the  variants  for  which 
OPARS  is  able  to  provide  flight  plans.  The  selected  aircraft 
also  comprise  over  70  percent  of  the  legs  computed  by  OPARS. 
The  significance  of  the  chosen  aircraft  is  also  apparent  when 
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considering  the  savings  that  FNOC  has  determined  from  using 
OPARS.  The  selected  aircraft  account  for  over  80  percent  of 
the  fuel  savings  as  calculated  by  FNOC.  The  FNOC  formula  used 
in  calculating  this  savings  will  be  discussed  later.  Table  IV 
shows  the  aircraft  that  are  used  in  this  analysis. 


Table  IV 

ANALYSIS 


AIRCRAFT   USED   AS    BASIS    FOR 


T4  3 

P3C 

C9D 

C20D 

UC12 

C9B 

P3A 

P3B 

HC130 

DC9 

KC130 

A  means  of  determining  the  route  that  the  unaided 
flight  planner  could  reasonably  be  expected  to  choose  in  lieu 
of  the  optimum  route  must  be  determined.  If  it  is  valid  to 
conclude  that  the  flight  planner  would  not  be  far  off  in  his 
estimate  of  the  optimum  altitude  at  which  to  fly,  then  an 
altitude  of  4  000  feet  on  either  side  of  the  optimum  should 
include  even  the  most  uneconomical  of  plans  that  the  flight 
planner  would  choose.  This  is  reasonable  due  to  extensive 
training  that  pilots  receive  in  flight  planning  and  their 
intimate  knowledge  of  their  aircraft.  It  has  been  indicated, 
that  for  the  S-3  aircraft  [Reference  6]  pilots  typically 
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select  an  altitude  as  much  as  10,000  feet  below  the  optimum 
altitude.  Since  information  of  this  type  is  not  available  for 
all  aircraft,  the  value  of  4000  feet  was  chosen  here.  This 
will  tend  to  under  estimate  the  savings  by  OPARS  if  indeed  the 
pilot's  range  of  error  is  greater  than  4000  feet.  For  the 
purposes  of  this  analysis  it  will  be  assumed  that  the  manual 
flight  planner  will  choose  a  flight  plan  that  is  the  optimum 
route  for  the  altitude  chosen.  The  overall  route,  as  far  as 
way  points  chosen,  will  be  the  same,  but  the  route  between  way 
points  will  be  allowed  to  vary  in  order  to  optimize  the  route 
at  each  altitude.  OPARS  will  select  the  optimum  jet  route 
between  user  way  points.  These  jet  routes  may  differ  at 
different  altitudes. 

Finally,  OPARS  is  capable  of  calculating  the  fuel 
required  to  fly  the  alternate  flight  plans  that  could  be 
chosen  by  the  flight  planner.  This  is  necessary  for  the 
comparisons  to  be  conducted  in  the  model  that  will  be 
discussed  later.  The  fuel  that  OPARS  calculated  to  be  used  on 
the  alternate  routes  will  be  optimum  for  that  route.  It  is 
unlikely  that  the  flight  planner  would  be  this  accurate  in 
his/her  calculations.  Therefore,  this  will  be  a  lower  bound 
on  the  percentage  of  additional  fuel  that  the  manual  planner 
would  require  for  the  flight. 

Information  is  not  available  on  the  rate  of  acceptance 
of  OPARS  flight  plan  recommendations.  Independent  studies 
from  two  aircraft  communities  [References  6  and  7]  have  shown 
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faith  in  the  system.  Additionally,  all  communities  have 
submitted  estimates  of  the  fuel  savings  that  they  feel  are 
afforded  by  use  of  OPARS .  Their  estimates  are  used  by  FNOC  in 
its  calculation  of  the  benefit  of  OPARS  from  fuel  savings.  It 
will  be  assumed  for  this  analysis  that  all  flights  are  flown 
using  an  OPARS  recommendation.  The  approach  used  in  the 
analysis  of  the  OPARS  system  is  to  determine  what  could  be 
saved  by  the  model.  Because  of  the  short  duration  of  the 
flights  and  the  ability  to  obtain  a  routing  just  minutes  prior 
to  the  actual  flight  time,  it  can  be  assumed  that  perfect 
information  is  available.  The  main  factor  that  would  decrease 
the  realized  savings  is  the  accuracy  of  the  weather 
information  itself.  For  OPARS,  it  will  be  assumed  that  the 
weather  model  used  by  FNOC  is  accurate. 

Currently,   FNOC   uses   the   following   formula   to 
calculate  the  savings  from  OPARS: 

Lx0.7Afx6=F 

where : 

•  L  is  the  number  of  legs  flown, 

•  M  is  the  maximum  internal  fuel  load, 

•  0  is  the  percent  of  fuel  estimated  to  be  saved  by  OPARS, 
and 

•  F  is  the  total  fuel  saved,  in  pounds. 

This  formula  has  three  major  faults.   The  first  is 
that  9  is  an  estimate  provided  by  the  squadrons  and  has  no 
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underlying  analysis.  Secondly,  not  all  flight  plans  are 
operational  plans,  as  is  assumed  by  the  above  equation.  A 
portion  of  the  plans  submitted  to  OPARS  are  duplicates  or  are 
for  more  than  72  hours  in  the  future.  The  third  fault  is  that 
the  formula  assumes  that  all  legs  are  loaded  to  70  percent  of 
maximum  internal  load.  With  increased  pilot  awareness  of  fuel 
conservation,  it  is  felt  that  this  is  too  high.  It  will  be 
shown  that  a  lower  figure  should  be  used. 

The  proposed  solution  to  this  formula  is, 

Lxi|rxMx6*xP=F 

where  in  addition  to  the  FNOC  formula: 

•  \j/  is  the  new  value  for  percent  fuel  load, 

•  0*  is  the  new  estimate  of  savings  and, 

•  P  is  the  probability  that  the  plan  is  an  operational  plan. 

Without  OPARS,  flight  plans  would  have  to  be  manually 
planned.  In  order  to  determine  what  the  unaided  flight 
planner  would  choose  for  a  flight  plan,  a  range  of  altitudes 
must  be  decided  upon.  As  earlier  described,  a  value  of  4000 
feet  on  either  side  of  the  optimum  was  chosen. 

Actual  flight  plan  requests  were  captured  for  a  forty- 
eight  hour  period.  These  plans  were  then  modified  to  force 
OPARS  to  compute  the  fuel  necessary  to  fly  at  specific 
altitudes.  In  this  way,  the  amount  of  fuel,  over  the  optimum, 
necessary  to  fly  at  the  various  altitudes  can  be  computed. 
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Once  the  minimum  amount  over  optimum  is  determined,  a  mean  and 
standard  deviation  from  all  aircraft  types  can  be  found. 

The  choices  that  the  unaided  flight  planner  will  make 
will  be  normally  distributed  about  this  mean  out  to  4000  feet 
on  either  side.  The  normal  distribution  function  combined 
with  the  fuel  use  curve  to  be  developed  and  the  number  of 
flight  plans  generated  will  provide  an  estimate  of  the  savings 
from  OPARS  over  manually  generated  flight  plans. 
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III.   RESULTS 

A.   OPTIMUM  TRACK  SHIP  ROUTING  SYSTEM 

To  determine  the  benefit  from  OTSR,  .the  shortest  path 
algorithm  and  the  greedy  algorithm  described  earlier  were 
used.  The  probability  of  wave  heights  greater  than  twelve 
feet  were  obtained  for  the  month  of  January  from 
climatological  charts  for  the  North  Atlantic  [Reference  3] . 
The  month  of  January  was  chosen  because  this  provided  the 
worst  case  for  sea  conditions.  This  would  provide  an  upper 
bound  on  the  damage  avoidance  estimate  for  OTSR. 

For  each  algorithm,  optimal  and  non-optimal,  a  route  from 
node  14  (075°W,  35°N)  to  node  157  (005°N,  50°W)  was  used. 
Once  the  difference  between  the  mean  values  for  the  damage 
sustained  on  the  optimal  and  non-optimal  routes  is  determined, 
it  will  be  multiplied  by  the  mean  number  of  routings  conducted 
per  month.  It  is  not  possible  to  accurately  determine  a  fuel 
savings  from  comparison  of  the  fuel  used  on  each  of  the 
routing  techniques  due  to  the  general  nature  of  the  fuel 
calculations.  By  using  the  fuel  required  by  the  general 
calculation,  it  was  possible  to  determine  the  most  economical 
path  under  each  of  the  routing  schemes  used.  If  the  only 
criterion  for  determining  the  route  had  been  damage  cost,  the 
algorithms  would  have  chosen  a  path  to  avoid  damage  even  at 
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the  expense  of  much  higher  fuel  consumption.  Clearly  this 
would  not  provide  the  optimum  route.  As  discussed  earlier,  it 
has  been  shown  by  other  studies  that  the  savings  in  fuel  from 
optimum  routing  of  ships  is  far  outweighed  by  that  of  damage 
avoidance . 

For  the  North  Atlantic,  99  routings  are  processed  per 
month.  As  discussed  earlier,  five  percent  (5)  of  these 
vessels  will  not  follow  the  OTSR  recommended  track,  and  of 
those  not  following  OTSR,  fifty  percent  (2.5)  will  encounter 
heavier  seas  and  possibly  sustain  damage.  Another  eleven 
percent  (10)  will  receive  rerouting  instructions  and  of  these, 
eighteen  percent  (2)  will  encounter  equal  or  heavier  seas.  In 
all,  4.5  ships  per  month  will  not  be  helped  by  the  OTSR 
system.  Therefore  the  estimate  of  the  savings  of  the  optimal 
over  non-optimal  routes  should  be  multiplied  by  94.5  vice  99, 
since  the  vessels  not  helped  by  OTSR  cannot  be  counted  as  a 
benefit  to  the  OTSR  system. 

A  similar  procedure  can  be  used  to  calculate  the  actual 
number  of  vessels  aided  by  OTSR  in  the  North  Pacific.  In  this 
case  the  number  of  vessels  should  be  126  vice  the  132  actually 
routed  by  the  Pearl  Harbor  Center.  Additionally,  since  the 
length  of  the  routes  are  typically  2.5  times  longer  in  the 
North  Pacific  than  in  the  North  Atlantic,  the  cost  of  the 
climatological  route  found  for  the  North  Atlantic  will  be 
multiplied  by  2.5  to  estimate  the  cost  of  a  climatological 
route  in  the  North  Pacific.   Once  a  yearly  savings  in  damage 
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avoidance  is  determined,  the  difference  between  the  damage 
known  to  occur  under  OTSR  and  that  incurred  under  optimal 
routing  in  known  weather  conditions  can  be  determined.  The 
results  of  this  analysis  using  the  route  given  above  is  shown 
in  Table  V. 


Table  V 

ROUTES . 


COSTS  ASSOCIATED  WITH  OPTIMAL  AND  NON-OPTIMAL 


Optimal 
Atl/Pac 

Non-Opt 
Atl/Pac 

Damage  Cost /Voyage 

433/1082 

2166/5415 

Ship  Days  Lost 

30.6 

152.7 

Number  of  Incidents/Year 

49.7 

248.4 

Damage  Cost /Year   Total 

2,127,006 

10,643,724 

Cost  of  Ship  Days  Lost/Year 

958,147 

4,781,342 

Total  Cost 

3,085,153 

15,425,066 

The  number  of  incidents  per  year  in  Table  V  was  determined 
using  the  rate  of  incidents  from  the  two  algorithms:  0.01  for 
the  optimal,  and  0.05  for  the  non-optimal.  Ship  days  lost 
were  determined  using  the  ratio  ship  days  lost  to  number  of 
incidents  from  Table  II.  The  cost  of  ship  days  lost  per  year 
is  computed  using  the  cost  of  MSC  vessels  per  day  at  31,312 
dollars  [Reference  8] . 

As  shown  in  Table  II  and  Table  V,  the  cost  of  damage  to 
ships  with  OTSR  recommendations  and  the  cost  of  damage  on 
optimally  routed  ships  is  very  close.  Once  the  cost  of  damage 
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from  Table  II  is  deducted  from  the  total  cost  of  non-optimal 
routing,  since  this  amount  of  damage  will  occur  with  or 
without  OTSR,  what  remains  is  the  savings  attributed  to  OTSR. 
That  savings  is  11,684,066  dollars. 

B.   OPTIMUM  PATH  AIRCRAFT  ROUTING  SYSTEM 
1.   Discussion  of  Procedure 

To  determine  the  savings  in  fuel  costs  afforded  by 
OPARS,  it  is  necessary  to  determine  the  cost  of  non-optimal 
routes.  As  discussed  earlier,  this  is  accomplished  using 
modified  flight  plans  and  resubmitting  them  to  OPARS. 
Original  flight  plans  were  obtained  as  they  were  submitted  by 
users  to  the  OPARS  model.  In  all,  364  flight  plans  were 
collected.  After  review,  it  was  determined  that  223  of  these 
were  in  fact  unique  flight  plans.  The  remainder,  upon  close 
examination,  were  found  to  be  in  one  of  the  following 
categories : 

•  Duplicates, 

•  Slight  modifications  of  a  basic  plan,  or 

•  Requested  for  more  than  72  hours  in  the  future. 

Duplicate  plans  could  be  readily  eliminated.  It  is 
unknown  why  they  were  submitted,  but  it  is  assumed  to  be  due 
to  user  impatience:  at  times,  the  queue  of  flight  plans 
submitted  may  become  long,  and  the  user  may  feel  that  his 
flight  plan  was  not  properly  submitted,  so  he  resubmits  it. 
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Flight  plans  that  are  slight  modifications  of  another 
plan  are  much  more  difficult  to  eliminate.  It  is  felt  that 
these  plans  are  an  attempt  by  the  user  to  experiment  and 
provide  a  range  of  options  available  for  the  requested  flight. 

The  third  category  was  eliminated  because  it  is  felt 
that  these  plans  would  be  resubmitted  at  a  future  date  that  is 
closer  to  the  actual  flight  time.  In  this  way,  actual  rather 
than  climatological  weather  would  be  used.  Flight  plans 
submitted  for  a  flight  time  greater  than  72  hours  in  the 
future  can  only  be  used  for  planning  purposes  and  cannot  be 
considered  as  operational  plans. 

A  summary  of  these  results  is  shown  in  Figure  3 .  The 
labeled  aircraft  are  those  that  were  used  in  the  analysis  of 
OPARS. 

Of  the  flight  plans  considered  operational,  57  where 
chosen  at  random  from  the  eleven  aircraft  types  chosen  for  the 
study.  These  flight  plans  were  then  modified  to  force  OPARS 
to  calculate  the  fuel  required  to  fly  at  each  of  four 
different  altitudes  evenly  spaced  over  the  range  of  altitudes 
available  to  the  aircraft.  In  this  way,  the  fuel  required  for 
non-operational  flight  plans  could  be  determined. 

The  altitude  restrictions  that  were  used  are  also 
available  to  the  users.  Twenty-one  of  the  original  flight 
plans  contained  constraints  on  the  altitude,  either  as  an 
upper  altitude  or  lower  altitude  restriction.  These 
constraints  were  removed  and  the  flight  plans  were  resubmitted 
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Figure  3     Operational  and  Non-operational  Flight 
Plans  from  the  Sample  Taken  of  User  Inputs 

along  with  the  original  and  modified  plans  to  provide  an 
alternate  base  line  for  the  optimum.  In  all,  306  flight  plans 
were  resubmitted  to  OPARS .  Because  of  weight  and  climb 
restrictions,  not  all  flight  plans  could  be  processed.  Of 
those  submitted  184  were  successfully  processed  and  provided 
355  individual  flight  legs  for  analysis.  Each  flight  plan 
could  contain  up  to  six  legs . 

When  OPARS  provides  the  completed  flight  plan,  three 
alternate  altitudes  and  the  fuel  required  for  that  altitude 
are  also  provided  for  each  leg  of  the  flight.  These  non- 
optimal  fuel  requirements  were  combined  with  the  fuel 
requirements  from  the  modified  flight  plans. 
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The  optimum  fuel  required  for  each  leg  was  compared  to 
all  available  non-optimal  fuel  requirements  for  that  same  leg. 
A  percentage  of  additional  fuel  required  for  each  altitude  was 
determined.  All  altitudes  and  the  percentage  of  additional 
fuel  required  for  that  altitude  within  4000  feet  of  the 
optimum  altitude  were  retained.  The  altitude  values  were  then 
coded  to  their  distance  from  the  optimum  and  are  displayed  in 
Figure  4.  By  coding  the  data,  it  was  possible  to  compare 
different  legs  for  different  aircraft  and  altitudes  on  a 
common  ground . 

As  expected,  the  additional  fuel  required  for  each  leg 
increases  as  the  distance  from  the  optimal  altitude  increases. 
Also  apparent  from  Figure  4  is  that  the  minimum  lies  to  the 
right  of,  i.e.,  at  a  higher  altitude  than,  the  optimum 
computed  by  OPARS .  On  further  inspection,  it  was  found  that 
this  was  due  to  the  altitude  restrictions  imposed  by  the  user. 
These  constraints  inhibited  OPARS  from  selecting  the  optimum 
altitude.  Figure  5  shows  the  same  information  as  Figure  4, 
but  in  this  case,  the  fuel  comparisons  were  made  against  the 
unconstrained  flight  plans  vice  the  original  as  entered  by  the 
user. 

The  minimum  to  the  right  of  the  optimum  still  persists 
in  Figure  5,  but  to  a  lesser  extent.  This  appears  to  be  due 
to  OPARS  reluctance  to  change  altitude  for  a  short  leg  if  the 
preceding  and  succeeding  legs  are  at  the  same  altitude,  so  as 
to  keep  the  flight  plan  at  a  level  altitude. 
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Figure  4     Percentage  of  Fuel  Over  Optimum  as 
Compared  to  Original  Flight  Plans 
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Figure  5     Percentage  of  Fuel  Over  Optimum  as 
compared  to  Unconstrained  Flight  Plans 
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For  the  purpose  of  flight  planning  and  navigation, 
altitude  is  not  continuous,  but  discrete  in  thousand  foot 
increments.  For  example  an  aircraft  is  assumed  to  be  at  15000 
feet  if  his  actual  altitude  is  between  14500  and  15500  feet. 
It  is  assumed,  as  discussed  earlier,  that  the  unaided  flight 
planner  would  choose  an  altitude  that  is  within  4000  feet  of 
the  optimum  altitude.  It  is  assumed  that  these  choices  will 
be  normally  distributed  about  the  optimum  with  a  mean  of  zero 
and  a  standard  deviation  of  4000.  By  converting  the  distance 
from  the  optimum  altitude  to  a  standard  normal,  and  using  the 
standard  normal  distribution  function,  the  probability  that 
the  unaided  flight  planner  will  choose  a  specific  incremental 
altitude  can  be  determined.  The  probability  of  choosing  a 
specific  altitude  and  the  mean  percent  fuel  required  above 
optimum  at  that  altitude  are  shown  in  Table  VI. 

Case  1  compares  the  modified  flight  plans  to  the 
original  flight  plans.  Case  2  compares  the  modified  flight 
plans  to  original  flight  plans  but,  with  the  altitude 
constraints  removed. 

The  expected  savings  in  terms  of  percentage  of  fuel 
used  over  the  optimum  will  be, 

SPax0a,   a=-4,-3, . . .,3,4. 

where , 

•  a  is  the  altitude  in  thousands  of  feet  from  the  optimum, 


35 


•  Pa  is  the  probability  that  a  specific  altitude  is  chosen, 
and 

•  0a  is  the  percentage  of  fuel  used  over  optimum  at  that 
altitude. 


Table  VI 


OPARS  SAVINGS  OVER  UNAIDED  FLIGHT  PLANS 


Alt,  x 
from 
Opt. 
ft  xlOO 

Mean  Percent 

Fuel  Used  Over 

Optimum 

Prob. 
L<x<U 
L= Lower 
U=Upper 

Percent  fuel 

Over  Optimum 

Given 

Altitude 

Case  1 

Case  2 

Case  1 

Case  2 

x<-35 

5.169 

5.294 

0.1908 

0.986 

1.010 

-35<x 
<-25 

3.022 

2.366 

0.0752 

0.227 

0.178 

-25<x 

<-15 

3.81 

4.068 

0.0879 

0.335 

0.358 

-15<x 
<-5 

1.547 

1.746 

0.0964 

0.149 

0.168 

-5<x<5 

0.738 

1.060 

0.0995 

0.073 

0.105 

5<x<15 

-0.950 

-0.022 

0.0964 

-0.092 

-0.002 

15<x<25 

1.156 

3.402 

0.0879 

0.102 

0.299 

2  5<x<35 

-0.400 

3.832 

0.0752 

-0.003 

0.288 

x>3  5 

2.608 

9.545 

0.1907 

0.497 

1.820 

Total 

2.324 

4.224 
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The  value  for  0*  in  the  following  equation, 

Lxi|rxATx0*xP=F/ 


where 


•  \|f  is  the  new  value  for  percent  fuel  load, 

•  6*  is  the  new  estimate  of  savings  and, 

•  P  is  the  probability  that  the  plan  is  an  operational  plan, 

as  computed  from  the  previous  equation,  is  then  4.224  for  case 
2  or  2.324  for  case  1  as  an  alternate.  The  value  of  4.224  is 
what  would  be  saved  by  OPARS  if  the  system  were  used  without 
altitude  constraints  allowing  OPARS  to  choose  the  optimum 
altitude  without  operator  intervention. 

As  can  be  seen  in  Table  VI,  nearly  67  percent  of  the 
savings  estimate  comes  from  the  tails  of  the  altitude 
distribution.  At  each  end  of  the  4000  foot  range,  the 
cumulative  probability  remaining  in  the  tails  is  great,  as  is 
the  percentage  of  fuel  used  over  the  optimum.  As  stated 
earlier,  the  choice  of  4000  feet  on  either  side  of  the  optimum 
ensured  that  from  available  data,  the  true  estimate  of  savings 
would  be  greater  than  the  value  determined  here. 

Currently,  when  computing  OPARS  fuel  savings,  FNOC 
uses  0.7  0  for  *F,  the  percentage  of  maximum  fuel  load.  Prior 
to  the  introduction  of  OPARS,  it  was  routine  to  load  aircraft 
to  100  percent  of  internal  fuel  load  for  every  flight.  When 
OPARS  was  introduced,  the  value  of  0.7  0  was  chosen  to  reflect 
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increased  awareness  of  fuel  conservation  and  to  underestimate 
the  savings  of  OPARS .  With  the  further  increases  in  fuel 
conservation  by  squadrons  and  the  need  to  maximize  the  number 
of  hours  flown  with  the  fuel  at  hand,  it  is  felt  that  this 
value  should  be  lowered  to  0.40.  To  explain  this  further,  it 
must  be  understood  that  in  order  to  determine  the  savings  of 
OPARS,  it  is  necessary  to  analyze  the  savings  at  the  level  of 
the  flight  leg.  The  current  method  assumes  that  each  leg  is 
loaded  to  7  0  percent  of  maximum  internal  load.  For  this  to  be 
the  case,  refueling  would  have  to  be  conducted  on  each  leg  of 
the  flight  plan.  This  is  not  so.  The  mean  loading  by  OPARS 
on  an  individual  leg  is  25  percent  of  internal  capacity  even 
taking  refueling  into  account.  Given  today's  concerns  over 
fuel  usage  and  conservation,  the  40  percent  chosen  here  for  *F 
is  felt  to  be  an  accurate  estimate  of  fuel  loading  in  the 
absence  Of  OPARS. 

Finally  the  probability,  P,  that  a  flight  plan  is  an 
operational  plan  is  determined  from  the  number  of  operational 
plans  observed  in  the  sample.  In  the  sample,  61  percent  of 
the  plans  were  operational.  For  this  analysis,  P  will  be  set 
at  0.75  to  eliminate  the  possibility  of  undercounting  the 
number  of  operational  plans. 


2 .   The  Savings  in  Fuel  from  the  use  of  OPARS 

For  this  analysis,  the  mean  number  of  legs  flown  each 
month  by  aircraft  type  in  1991  was  used  to  arrive  at  the 
savings  in  fuel  by  OPARS. 

Using  the  revised  estimation  procedure  for  each 
aircraft  type,  OPARS  is  estimated  to  save  6.773  million 
dollars  when  using  the  preferred  case  2  data,  and  3.726 
million  dollars  if  the  case  1  data  is  used.  This  is  in 
comparison  to  FNOC's  estimate  of  8.348  million  dollars. 
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IV.   SUMMARY  AND  RECOMMENDATIONS 

A.   OPTIMUM  TRACK  SHIP  ROUTING  SYSTEM 

Calculations  to  determine  the  savings  realized  by  the  U. 
S.  Navy  from  OTSR  are  based  on  two  related  algorithms.  The 
first  uses  an  adaptation  of  Dijkstra's  algorithm  to  determine 
the  shortest  path  across  the  North  Atlantic  with  all  wave 
heights  known.  The  second  algorithm  uses  a  greedy  strategy 
and  looks  only  at  the  next  accessible  nodes  that  can  lead  to 
the  desired  destination  and  chooses  the  least  cost  of  those 
available. 

Surprisingly,  it  was  found  that  the  annual  damage  costs 
sustained  under  the  first  algorithm  closely  match  the  damage 
costs  that  are  experienced  under  the  OTSR  system.  The 
significance  of  this  is  not  explored  here,  but  it  may  be 
possible  to  show  that  the  OTSR  route  is  quite  close  to  the 
true  optimum  route.  The  second  algorithm  was  used  to 
determine  damage  costs  in  the  absence  of  OTSR.  In  this  case 
Commanding  Officers  would  be  required  to  rely  on 
climatological  or  short  range  forecasts  to  choose  their  route. 
As  has  been  shown,  this  results  in  much  greater  damage  costs. 

To  determine  an  estimate  of  the  savings  from  OTSR,  a 
single  route  from  the  northeastern  coast  of  the  United  States 
to  the  southern  tip  of  England  was  used.   Random  numbers 
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chosen  from  a  uniform  distribution  determined  wave  heights  and 
probability  of  damage  at  each  node.  The  same  wave  heights  and 
probabilities  were  used  for  each  of  the  algorithms  by  using 
the  same  seed.  Each  algorithm  was  run  one  hundred  times  with 
different  sea  conditions  and  damage  probabilities  to  determine 
a  mean  cost  for  the  route  under  optimal  and  non-optimal 
routing. 

This  mean  of  the  non-optimal  routes  was  then  multiplied  by 
the  number  of  routings  per  year  in  the  North  Atlantic  and  the 
North  Pacific.  In  the  North  Pacific  an  additional  scaling 
factor  was  used  due  to  the  length  of  voyages  there.  The 
annual  cost  of  non-optimal  routing  once  decreased  for  damage 
costs  occurring  even  while  under  OTSR  control  was  11.7  million 
dollars . 

B.   OPTIMUM  PATH  AIRCRAFT  ROUTING  SYSTEM 

This,  the  second  product  under  study,  was  evaluated  using 
the  OPARS  model  itself.  Modified  flight  plans  were 
resubmitted  to  OPARS  to  determine  the  amount  of  fuel  required 
for  a  non-optimal  flight.  The  flight  plans  had  originally 
been  copied  as  the  requests  were  received.  They  were  then 
modified  to  require  OPARS  to  determine  fuel  loading  if  the 
plan  were  flown  at  a  specific  altitude. 

Eleven  aircraft  types  were  chosen  and  numerous  flight 
plans  from  each  were  modified.  Four  altitudes  were  chosen  for 
each  aircraft,  depending  on  its  capabilities,  in  order  to 


41 


bracket  the  optimum  altitude  computed  for  the  original.  For 
each  flight  plan  then,  five  flight  plans  were  resubmitted;  the 
original  and  four  at  modified  altitudes. 

Once  the  fuel  required  for  each  altitude  was  determined, 
it  was  compared  to  the  optimum.  In  this  way,  a  percentage  of 
fuel  required  over  optimum  could  be  computed.  The  distance 
from  the  optimum  altitude  was  also  determined.  The  percentage 
of  fuel  over  optimum  was  plotted  against  the  distance  from  the 
optimum  altitude. 

The  percentage  of  fuel  over  optimum  was  multiplied  by  the 
probability  of  a  manual  flight  planner  choosing  that  altitude. 
The  probability  of  choosing  a  particular  altitude  was  based  on 
a  normally  distributed  random  variable  with  mean  zero  and 
standard  deviation  4000  feet.  In  doing  this,  an  aggregate 
value  for  the  estimate  of  fuel  saved  by  OPARS  was  found. 

Once"  this  value  was  entered  into  the  modified  FNOC  fuel 
equation,  a  fuel  savings  of  6.8  million  dollars  was  estimated. 
The  modified  FNOC  equation  changes  several  of  the  parameters 
used  by  FNOC  in  their  current  calculation.  First,  the  savings 
estimate  described  above  is  used  instead  of  an  estimate 
provided  by  each  aircraft  squadron.  Next  a  new  value  for  the 
percent  of  maximum  fuel  load  that  would  be  loaded  in  the 
absence  of  OPARS  is  used.  Currently  this  value  is  at  7  0 
percent.  Actual  loading  by  OPARS  is  25  percent  based  on  the 
observed  flight  plans.  A  value  of  40  percent  was  used  in  the 
modified  equation.    Finally  a  parameter  to  indicate  the 
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probability  that  the  flight  plan  was  actually  an  operational 
plan  was  added.  A  value  of  0.75  was  chosen  for  this  parameter 
based  on  the  observed  flight  plans. 

C.   RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

As  has  been  shown,  FNOC  does  indeed  provide  a  valuable 
service.  The  cost  savings  demonstrated  here  involved  only  two 
of  the  many  products  produced  by  FNOC.  Further  study  should 
be  devoted  to  quantifying  the  remaining  products  not  covered 
here.  Work  must  also  be  done  to  accurately  determine  the  fuel 
savings  that  can  be  attributed  to  OTSR. 

Additionally,  work  is  being  done  to  build  a  computer  model 
that  could  be  placed  aboard  ships  to  generate  their  own 
optimum  routes.  When  this  is  accomplished,  studies  should  be 
done  to  determine  the  added  benefit  from  having  this 
capability  aboard  ships. 

In  order  to  obtain  more  exact  estimates  of  the  savings 
from  damage  avoidance  provided  by  OTSR,  follow  up  to  damage 
reports  submitted  to  the  Naval  Safety  Center  must  be 
conducted.  The  estimates  provided  here  appear  to  be  low  and 
should  be  revised  to  obtain  a  more  accurate  benefit  from  OTSR. 

The  procedures  used  here  for  both  OTSR  and  OPARS  provide 
the  framework  for  further  study.  In  each  case  an  estimation 
of  the  savings  provided  by  the  product  is  given.  Further 
study  should  be  given  to  sensitivity  analysis  of  the 
parameters  involved.    Namely,   in  the  case  of  OPARS,  the 
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percentage  of  fuel  used  over  optimum  and  the  percentage  of 
maximum  internal  load  should  be  studied.  Additionally,  it 
will  be  necessary  to  more  accurately  determine  the  range  of 
altitudes  that  would  be  chosen  by  manual  flight  planners.  For 
OTSR,  a  more  accurate  method  of  calculating  fuel  use  should  be 
investigated  to  determine  more  accurately  any  benefit  gained 
from  fuel  savings  in  optimum  routing. 

Additionally,  in  the  case  of  OTSR,  actual  routings  and 
weather  conditions  should  be  collected  for  use  with  the 
modified  Dijkstra  algorithm  to  determine  possible  future  gains 
for  OTSR  savings  should  forecasting  methods  improve. 
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APPENDIX  A 

Optimal  Path  Algorithm 

This  program  reads  in  a  sparse  graph  and  determines  the 
shortest  path  through  the  graph. 


CONST  MAX=160; 
START=14; 
STOP=157; 
RANDSEED=7654321; 

TYPE  PADJ="ADJNODES; 
PEDGE= "EDGES; 
EDGES =RECORD 

REC : PAD J ; 
NEXTNODE : PEDGE ; 
PARENT : PAD J ; 
WT : REAL ; 
END; 
AD JNODES  =RECORD 

NODE : INTEGER ; 
NEXTNODE : PEDGE ; 
NODE} 

DIST:REAL; 
PWAVE : REAL ; 
PDAM:REAL; 
PRED:PADJ; 
DAM : REAL ; 
END; 
HE ADNODE= ARRAY [1 . .MAX] 
QUEUE = ARRAY [ 1 .  . MAX ]  OF  PAD J ; 
AY= ARRAY [ 1 . . MAX ]  OF  REAL ; 
7]  OF  REAL; 


{POINTS  TO  THE  TAIL} 
{NEXT  ADJACENT  ARC} 
{POINTS  TO  THE  HEAD} 
{WEIGHT  OF  ARC} 


{NODE  NUMBER} 

{POINTS  TO  NEXT  ADJACENT 

{DISTANCE  FROM  SOURCE} 

{ PROB  OF  A  WAVE  HT . } 

{ PROB  DAMAGE  GIVEN  WAVE  HT . } 

{POINTER  TO  PREDECESSOR} 


OF  PADJ; 


PC=ARRAY[1 
G=RECORD 


GRAPH : HEADNODE ; 
LENGTH: INTEGER ; 


{THE  GRAPH} 

{THE  LENGTH  OF  THE  GRAPH} 


END; 
Q=RECORD 


END; 


PRIQ : QUEUE ; 
SIZE: INTEGER 


VAR  I, J: INTEGER; 
DATAOUT : TEXT ; 
GRAPH1:G; 
PI , P2 , PNTLAT , PNTLON : AY ; 

Q1:Q; 

PDAMAGE , C DAMAGE : PC ; 
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CPNT:PADJ; 
CP:PEDGE; 
TOTAL : REAL ; 

PROCEDURE  RPROB(VAR  FIRST: AY;VAR  SEC: AY); 

{READS  PROBABILITY  OF  WAVE  HEIGHT  FROM  INPUT  FILE} 

VAR  I : INTEGER; 
•   DATAIN : TEXT ; 

BEGIN 

ASSIGN (DATAIN, 'C:\PASCAL\WAVES.TXT' ) ; 

RESET (DATAIN) ; 

FOR  I:=l  TO  MAX  DO 

READLN (DATAIN, FIRST [I] ,SEC[I] ) ; 
END; 

PROCEDURE  RLATLON(VAR  FIRST: AY; VAR  SEC: AY); 
{READS  LAT  AND  LONG  FROM  INPUT  FILE} 

VAR  I : INTEGER ; 

DATAIN: TEXT; 

BEGIN 

ASSIGN (DATAIN, 'C:\PASCAL\LATLON.TXT' ) ; 

RESET (DATAIN) ; 

FOR  I:=l  TO  MAX  DO 

READLN (DATAIN, FIRST [I] ,SEC[I] ) ; 
END; 

FUNCTION" FUEL ( BEG : INTEGER ; EN : INTEGER ) : REAL ; 
{COMPUTES  FUEL  REQUIRED  FOR  AN  ARC} 

VAR  LATS , LATD , LONS , LOND , COSDI ST , J , TEMP : REAL ; 

BEGIN 

LATS :=PNTLAT [BEG]* (PI/180) ; 
LONS :=PNTLON [BEG] * (PI/180)  ; 
LATD : =PNTLAT [ EN] * ( PI / 1 8  0 )  ; 
LOND:=PNTLON[EN] * (PI/180) ; 
COSDIST:=( SIN (LATS) *SIN(LATD) )+ (COS (LATS) *COS(LATD) * 
(COS (LOND -LONS) ) ) ; 
J:=0; 

TEMP: =1.0; 
REPEAT 
BEGIN 

J:=J+0.01; 
TEMP:=COS(J) ; 
END; 
UNTIL ( (COSDIST>=TEMP)OR(J=3 . 14 ) ) ; 
FUEL:=( ( (J*180/PI)*60)/15)* (650/24) ; 
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END; 

FUNCTION  BUILDWTS (VAR  DEST: PEDGE) : REAL; 

{BUILD  THE  WEIGHT  OF  THE  ARC  FROM  FUEL  AND  DAMAGE} 

VAR  I : INTEGER ; 

COST, DIST: REAL; 

BEGIN 

COST:=0; 

IF( (DEST". REC ".PWAVE<=P2 [DEST" .REC" .NODE] ) AND 
(P2 [DEST". REC". NODE] <1) ) THEN 
FOR  I:=l  TO  3  DO 

IF (DEST" . REC" . PDAM< = PDAMAGE [ I ] ) THEN 
COST : =CDAMAGE [  I  ]  ; 
IF ( (COST=0 ) AND (DEST" .REC" . PDAM<=P1 [DEST" .REC" .NODE] ) AND 
(PI [DEST". REC". NODE] <1) ) THEN 
FOR  I: =4  TO  5  DO 

IF  (DEST".REC".PDAM<=PDAMAGE[I] ) THEN 
COST : =CDAMAGE [ I ] ; 
IF(COST=0)THEN 

FOR  I: =6  TO  7  DO 

IF  ( DEST". REC ".PDAM<= PDAMAGE [I] ) THEN 
COST : =  CDAMAGE [  I  ]  ; 
BUILDWTS : =COST+ (FUEL (DEST" . PARENT" . NODE , 
DEST". REC ".NODE) *45.86) ; 
END; 

PROCEDURE  HEAPIFY (VAR  NQ : Q; START : INTEGER) ; 

VAR  SMALLEST , L , R : INTEGER ; 
TEMP : PAD J; 

BEGIN 

L:=2*START; 

R:=(2*START)+1; 

IF  ((L<NQ.SIZE)  AND 

(NQ.PRIQ[L] ". DIST<NQ.PRIQ[ START] ".DIST) )  THEN 
SMALLEST :=L 
ELSE 

SMALLEST: =START; 
IF  ( (R<NQ.SIZE)   AND 

(NQ.PRIQ[R] ".DIST<NQ.PRIQ[ SMALLEST] ".DIST) )  THEN 
SMALLEST: =R; 
IF  (SMALLESToSTART)  THEN 
BEGIN 

TEMP  :  =NQ  .  PRIQ  [  START]  ; 
NQ  .  PRIQ  [  START]  :  =NQ  .  PRIQ  [  SMALLEST]  ; 
NQ . PRIQ [ SMALLEST] : =TEMP ; 
HEAPIFY (NQ, SMALLEST) ; 
END;    {IF  SWAPPED} 
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END;    {PROCEDURE  HEAPIFY} 

PROCEDURE  INSERTPQ ( PNTR : PADJ ;  VAR  NQ : Q ) ; 

VAR  I : INTEGER; 

BEGIN 

NQ . SIZE : =NQ . SIZE+1 ; 
I:=NQ.SIZE; 

WHILE  ((I>1)  AND  (NQ.PRIQ[I  DIV  2 ] " . DIST>PNTR" . DISTl 
DO 

BEGIN 

NQ . PRIQ [ I ] : =NQ . PRIQ [ I  DIV  2 ] ; 
I:=I  DIV  2; 
END;    {WHILE} 
NQ.PRIQ[I] :=PNTR; 
END;    {PROCEDURE  INSERTPQ} 

PROCEDURE  BUILDPQ(THISGRAPH:G;VAR  PQ:Q); 

VAR  I : INTEGER; 

BEGIN 

FOR  I:=l  TO  THISGRAPH. LENGTH  DO 
BEGIN 

INSERTPQ (THISGRAPH. GRAPH [I] ,  PQ)  ; 
END;    {FOR} 
END;    {PROCEDURE  BUILDPQ} 

FUNCTION  EXTRACTMIN (VAR  PQ : Q) : PADJ ; 

BEGIN       {FUNCTION  EXTRACTMIN} 

EXTRACTMIN : =PQ . PRIQ [ 1 ] ; 

PQ.PRIQ[1] :=PQ. PRIQ [PQ. SIZE] ; 

PQ.SIZE:=PQ.SIZE-1; 

HEAPIFY (PQ,1) ; 
END;    {FUNCTION  EXTRACTMIN} 

FUNCTION  EMPTYPQ(VAR  PQ:Q) : BOOLEAN; 

BEGIN 

IF  (PQ.SIZE=0)  THEN 

EMPTYPQ : =TRUE 
ELSE 

EMPTYPQ : = FALSE ; 
END;    {FUNCTION  EMPTYPQ} 

PROCEDURE  MAKEGRAPH(VAR  THISGRAPH :G) ; 
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{  THIS  PROCEDURE  GENERATES  THE  GRAPH.} 

TYPE  ARY=ARRAY [ 1 . . 2 ]  OF  INTEGER ; 

VAR  NEWREC : PEDGE ; 
CP, LP: PEDGE; 
DUPE , FOUND : BOOLEAN ; 
DATAIN : TEXT ; 
NEWNODE : ARY ; 
I, J, TEMP: INTEGER ; 
NEWWT: INTEGER; 

BEGIN 

ASSIGN (DATAIN, 'C:\PASCAL\G.TXT') ; 
RESET (DATAIN) ; 
NEWWT :=0; 

READLN ( DATAIN , THI SGRAPH . LENGTH ) ; 
WHILE (NOT  EOF (DATAIN))  DO 
BEGIN 

FOR  I:=l  TO  2  DO 

READ ( DATAIN, NEWNODE [ I] ) ; 
READLN (DATAIN) ; 

IF  (NEWNODE [1]<>NEWN0DE [2 ] )  THEN         {IF  NOT  A 
SELF  LOOP} 
BEGIN 

DUPE:=FALSE; 
FOUND :=FALSE; 

IF  ( THI  SGRAPH.  GRAPH  [NEWNODE  [1]  ]  * .  NEXTNODEoNIL  \ 
THEN  {FIRST  NODE} 
BEGIN 

LP:=NIL; 

CP : =THISGRAPH . GRAPH [NEWNODE [ 1 ] ] A .NEXTNODE ; 

REPEAT 

DUPE:  =  (CP/v.REC".NODE=NEWNODE[2]  )  ; 
FOUND : = ( CP^ . REC^ . NODE>NEWNODE [2 ] ) ; 
IF     (NOT (FOUND)    AND   NOT (DUPE))     THEN 
BEGIN 

LP:=CP; 

CP:=CPA. NEXTNODE;   • 
END; 
IF  (CP=NIL)  THEN 
FOUND : =TRUE 
UNTIL  (FOUND  OR  DUPE  OR  (CP=NIL) ) ; 
IF  (NOT(DUPE))  THEN 
BEGIN 

NEW  (NEWREC)  ; 

NEWREC" . NEXTNODE : =CP ; 

NEWREC". REC:=THI SGRAPH. GRAPH [NEWNODE [2]  ]  ; 

NEWREC" . PARENT : = 

THI SGRAPH . GRAPH [NEWNODE [ 1 ] ] ; 
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ARC} 


NEWREC ".WT:=0; 
IF  FOUND  THEN 

IF  (LP=NIL)  THEN 

THISGRAPH . GRAPH [NEWNODE [ 1 ] ] " . 
NEXTNODE : = NEWREC 
ELSE 

LP" . NEXTNODE : =NEWREC ; 
END;   {IF  NOT  DUPE} 
IF  DUPE  THEN         {ADDS  ONLY  SMALLEST 

IF  (NEWWT<CP/S.WT)  THEN     {MULTIPLE  ARCS} 
CP".WT:=NEWWT; 
END 
ELSE  {IF  FIRST  EDGE} 

BEGIN 

NEW(NEWREC) ; 

NEWREC ~ . NEXTNODE : =NIL ; 

NEWREC  A . REC : =THI SGRAPH . GRAPH [ NEWNODE [  2  ]  ]  ; 

NEWREC ~ . PARENT : =THI SGRAPH . GRAPH [NEWNODE [ 1] ] ; 

NEWREC  A.WT: =NEWWT ; 
THI SGRAPH . GRAPH [ NEWNODE [ 1 ] ] A . 
NEXTNODE : =NEWREC ; 
END; 
END;       {IF  NO  SELF  LOOP} 
END;       {WHILE} 
END;        {PROCEDURE  MAKEGRAPH} 


PROCEDURE  INITGRAPH(VAR  G1:G); 

BEGIN 

FOR  I:=l  TO  MAX  DO 
BEGIN 

NEW(G1.GRAPH[I] ) ; 
WITH  Gl. GRAPH [I]"  DO 
BEGIN 

NODE:=I; 
DIST:=10000000; 
PWAVE:=RANDOM; 
PDAM:=RANDOM; 
PRED:=NIL; 
NEXTNODE : =NIL ; 
DAM:=0.0; 
END; 
END; 
Gl. LENGTH :=0; 
END;      {PROCEDURE  INITGRAPH} 

PROC  EDURE  INI TQ ( VAR  TH I SQ : Q )  ; 
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VAR  I: INTEGER; 

BEGIN 

FOR  I :=1  TO  MAX  DO 
BEGIN 

NEW(THISQ.PRIQ[I]  )  ; 
THISQ.PRIQ[I] :=NIL; 
END; 
THISQ.SIZE:=0; 
END; 


PROCEDURE  DI JKSTRA ( VAR  THI SGRAPH : G ; SOURCE : INTEGER ) ; 

VAR  I : INTEGER; 
WT : REAL ; 
HERE : PEDGE ; 
THISNODE:PADJ; 

BEGIN 

THI SGRAPH . GRAPH [SOURCE] A . PRED : =THI SGRAPH . GRAPH [ SOURCE] ; 
THISGRAPH . GRAPH [ SOURCE ] A . DIST : =  0 ; 
BUILDPQ( THI SGRAPH, Ql) ; 
WHILE  (NOT(EMPTYPQ(Ql) ) )  DO 
BEGIN 

THISNODE:=EXTRACTMIN(Ql) ; 
HERE  :  =THISNODE/v .  NEXTNODE  ; 
WHILE  (HEREoNIL)  DO 
BEGIN 

WT : =BUILDWTS ( HERE ) ; 

IF(HERE".REC/s.DIST>THISNODE^.DIST+WT)  THEN 
BEGIN 

HERE^  .  REC"  .  DIST :  =THISNODE/N .  DIST+WT; 
HERE" .REC" . PRED : =THISNODE ; 
FOR  I:=( (Ql.SIZE+1)  DIV  2)  DOWNTO  1  DO 
HEAPIFY(Q1,I) ; 
END; 
HERE : =HERE" . NEXTNODE ; 
END; 
END; 
END;       {PROCEDURE  DI JKSTRA} 


BEGIN         {MAIN  PROGRAM} 

ASSIGN(DATAOUT, 'C:\PASCAL\0UTPUT3 .TXT' ) ; 

REWRITE (DATAOUT) ; 

RPROB ( PI , P2 ) ; 

RLATLON(PNTLAT, PNTLON) ; 

GRAPH 1. GRAPH [START] A.DIST:=0; 

PDAMAGE[1] : =0.858 ; 

PDAMAGE[2] :=0.21; 
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PDAMAGE [ 3 ] 
PDAMAGE [ 4 ] 
PDAMAGE [ 5 ] 
PDAMAGE [ 6 ] 
PDAMAGE [ 7 ] 
CDAMAGE [ 1 ] 
CDAMAGE [ 2 ] 
CDAMAGE [ 3 ] 
CDAMAGE [ 4 ] 
CDAMAGE [ 5 ] 
CDAMAGE [ 6 ] 
CDAMAGE [ 7 ] 
INITGRAPH ( 


=  0 
=  0 
=  0 
=  0 
=  0 


087; 
0354 
0027 
0024 
0003 
=340771 
=312196 
=312196 
=129969 
=129969 
=48427 
=48427 
GRAPH1) 


.PRED:=NIL; 
,DIST:=10000000 


MAKEGRAPH (GRAPH1 
INITQ(Ql) ; 
FOR  I:=l  TO  100  DO 
BEGIN 

FOR  J:=l  TO  MAX  DO 
BEGIN 

GRAPH 1. GRAPH [J] 
GRAPH 1. GRAPH [J] 
Q1.PRIQ[J] :=NIL; 
END; 
Q1.SIZE:=0; 

DIJKSTRA(GRAPH1, START) ; 
CPNT:=GRAPH1. GRAPH [STOP] ".PRED; 

TOTAL : =FUEL (CPNT". NODE , GRAPH 1 . GRAPH [STOP] ".NODE) *45.86 
WRITE (DATAOUT, STOP : 4 ) ; 
WHILE  (CPNT~  .  PRED"  .NODEoCPNT"  .NODE; 
BEGIN 
"  TOTAL := TOTAL + 

( FUEL ( CPNT" . PRED" . NODE , C  PNT" 
WRITE ( DATAOUT , CPNT" . NODE : 4 ) ; 
CPNT : =CPNT" . PRED ; 
END; 
WRITE ( DATAOUT , START : 4 ) ; 

WRITE (DATAOUT, GRAPH 1 .GRAPH [STOP] " .DIST: 10 : 2 ) ; 
WRITELN (DATAOUT, GRAPH 1 .GRAPH [STOP] * . DIST-TOTAL  : 10 : 2 ) 
FOR  J:=l  TO  MAX  DO 
BEGIN 

GRAPH 1. GRAPH [J] 
GRAPH 1. GRAPH [J] 
END; 
END; 

CLOSE (DATAOUT) ; 
END. 


DO 


NODE)*45.86) 


PWAVE := RANDOM; 
PDAM:=RANDOM; 
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APPENDIX  B 

Non-optimal  Routing  Algorithm 

This  program  reads  in  a  sparse  graph  and  determines  a 
non-optimal  path  through  the  graph  based  on  only  the  next 
immediately  available  nodes. 

CONST  MAX=160; 
START=14; 
STOP=157; 
RANDSEED=7 654321; 


TYPE  PADJ=~ADJNODES; 
PEDGE = "EDGES; 
EDGES=RECORD 

REC : PADJ ; 

NEXTNODE : PEDGE ; 

PARENT : PADJ ; 

WT : REAL ; 
END; 
AD JNODES  =RECORD 

NODE: INTEGER; 

NEXTNODE : PEDGE ; 

DIST:REAL; 

PWAVE : REAL ; 

PDAM : REAL ; 

PRED : PADJ ; 

DAM : REAL ; 

QIN : BOOLEAN; 
END; 
HE ADNODE = ARRAY [ 1 .  . MAX ]  OF 
AY= ARRAY [ 1 .  . MAX  ]  OF  REAL ; 
PC=ARRAY[1. .7]  OF  REAL; 
G=RECORD 

GRAPH : HEADNODE ; 

LENGTH : INTEGER ; 
END; 


{POINTS  TO  THE  TAIL} 
{NEXT  ADJACENT  ARC} 
{POINTS  TO  THE  HEAD} 
{WEIGHT  OF  ARC} 


{NODE  NUMBER} 
{POINTER  TO  NEXT  NODE} 
{DISTANCE  FROM  SOURCE} 
{PROB  OF  WAVE  HEIGHT} 
{PROB  OF  DAMAGE  GIVEN  HT.} 
{POINTER  TO  THE  PREDECESSOR} 
{DAMAGE  ENCOUNTERED} 
{IS  NODE  ON  PATH  TO  DEST.} 

PADJ; 


{THE  GRAPH} 

{THE  LENGTH  OF  THE  GRAPH} 


VAR  I , J: INTEGER ; 
DATAOUTrTEXT; 
GRAPH1 : G ; 

PI , P2 , PNTLAT , PNTLON : AY ; 
PDAMAGE , CDAMAGE : PC ; 
CPNT:PADJ; 
CP : PEDGE; 
TOTAL : REAL ; 

PROCEDURE  RPROB(VAR  FIRST: AY; VAR  SEC:AY); 

VAR  I : INTEGER; 
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DATAIN : TEXT ; 

BEGIN 

ASSIGN (DATAIN, 'C:\PASCAL\WAVES.TXT' ) ; 

RESET (DATAIN) ; 

FOR  I:=l  TO  MAX  DO 

READLN (DATAIN, FIRST [I] ,SEC[I] ) ; 
END; 

PROCEDURE  RLATLON(VAR  FIRST : AY; VAR  SEC: AY); 

VAR  I : INTEGER; 

DATAIN : TEXT ; 

BEGIN 

ASSIGN (DATAIN, 'C:\PASCAL\LATLON.TXT' ) ; 

RESET (DATAIN) ; 

FOR  I:=l  TO  MAX  DO 

READLN (DATAIN, FIRST [I] ,SEC[I] ) ; 
END; 

FUNCTION  FUEL (BEG: INTEGER; EN: INTEGER) :REAL; 

VAR  LATS , LATD , LONS , LOND , COSDI ST , J , TEMP : REAL ; 

BEGIN 

LATS :=PNTLAT[ BEG] * (PI/180)  ; 
LONS :=PNTLON[ BEG] * (PI/180) ; 
LATD:=PNTLAT[EN] * (PI/180) ; 
LOND:=PNTLON[EN] * (PI/180) ; 
COSDlST:=( SIN (LATS) *SIN(LATD) )+ 

(COS (LATS) *COS(LATD) * (COS (LOND -LONS) 
J:=0; 

TEMP : =1.0; 
REPEAT 
BEGIN 

J:=J+0.01; 
TEMP:=COS(J) ; 
END; 
UNTIL ( (COSDIST>=TEMP)OR(J=3 .14) ) ; 
FUEL:=( ( (J*180 /PI) *60) /15)* (650/24) ; 
END  ; 

FUNCTION  BUILDWTS(VAR  DEST : PEDGE) :REAL; 

VAR  M: INTEGER ; 

COST, DIST: REAL; 

BEGIN 

COST:=0; 

IF(  (DEST/S.REC/S.PWAVE<  =  P2  [DEST"  .REC"  .NODE]  ) 
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AND(P2  [DEST/V.REC/V.N0DE]<1)  )  THEN 
FOR  M:=l  TO  3  DO 

IF (DEST" .REC" . PDAM< = PDAMAGE [M] ) THEN 
COST : =CDAMAGE [M] ; 
IF(  (COST=0)AND(DEST/V„REC/S.PDAM<  =  P1  [DEST~ .REC" .NODE]  ) 
AND(P1[DEST".REC~.N0DE]<1) ) THEN 
FOR  M:=4  TO  5  DO 

IF  (DEST".REC".PDAM<=PDAMAGE[M] ) THEN 
COST : =CDAMAGE [M] ; 
IF(COST=0)THEN 

FOR  M:=6  TO  7  DO 

IF  (DEST".REC/S.PDAM<  =  PDAMAGE  [M]  )  THEN 
COST:=  CDAMAGE[M]; 
BUILDWTS : =COST+ 

( FUEL ( DEST" . PARENT" . NODE , DEST" . REC" . NODE ) *  4  5 . 8  6 )  ; 
END; 


PROCEDURE  MAKEGRAPH(VAR  THISGRAPH :G) ; 

{  THIS  PROCEDURE  GENERATES  THE  GRAPH.} 

TYPE  ARY=ARRAY [ 1 . . 2 ]  OF  INTEGER ; 

VAR  NEWREC : PEDGE ; 
CP, LP: PEDGE; 
DUPE , FOUND : BOOLEAN ; 
DATAIN : TEXT ; 
NEWNODE : ARY ; 
I, J, TEMP: INTEGER; 
NEWWT : INTEGER ; 

BEGIN 

ASSIGN (DATAIN, 'C:\PASCAL\G.TXT') ; 
RESET (DATAIN) ; 
NEWWT: =0 ; 

READLN( DATAIN, THISGRAPH. LENGTH) ; 
WHILE (NOT  EOF (DATAIN))  DO 
BEGIN 

FOR  I:=l  TO  2  DO 

READ ( DATAIN, NEWNODE [  I  ]  )  ; 
READLN (DATAIN) ; 
IF  (NEWNODE [1]<>NEWN0DE [2 ] )  THEN  {IF  NOT  A  SELF 


LOOP} 


BEGIN 

DUPE:=FALSE; 
FOUND :=FALSE ; 

IF  (THISGRAPH. GRAPH  [NEWNODE  [1]  ]  *  .NEXTNODEoNIL) 
THEN  {FIRST  NODE} 
BEGIN 

LP:=NIL; 
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CP:=THISGRAPH.GRAPH[NEWN0DE[1] ] " . NEXTNODE ; 
REPEAT 

DUPE:=(CP".REC".NODE=NEWNODE[2] ) ; 
FOUND:=(CP".REC".NODE>NEWNODE[2] ) ; 
IF  (NOT (FOUND)  AND  NOT (DUPE))  THEN 
BEGIN 

LP:=CP; 

CP : =CP" . NEXTNODE ; 
END; 
IF  (CP=NIL)  THEN 
FOUND : =TRUE 
UNTIL  (FOUND  OR  DUPE  OR  (CP=NIL) ) ; 
IF  (NOT (DUPE))  THEN 
BEGIN 

NEW(NEWREC) ; 

NEWREC" . NEXTNODE : =CP ; 

NEWREC/s.REC:=THISGRAPH.GRAPH[NEWNODE[2]  ]  ; 

NEWREC" . PARENT : = 

THISGRAPH . GRAPH [NEWNODE [ 1 ] ] ; 
NEWREC~.WT:=0; 
IF  FOUND  THEN 

IF  (LP=NIL)  THEN 

THISGRAPH . GRAPH [NEWNODE [ 1 ] ] A . NEXTNODE 

: =NEWREC 
ELSE 

LP" . NEXTNODE : =NEWREC ; 
END;   {IF  NOT  DUPE} 
IF  DUPE  THEN        {ADDS  ONLY  SMALLEST  ARC} 
IF  (NEWWT<CP/S.WT)  THEN    {MULTIPLE  ARCS} 
CP".WT:=NEWWT; 
END 
ELSE         {IF  FIRST  EDGE} 
BEGIN 

NEW(NEWREC) ; 

NEWREC" . NEXTNODE : =NIL ; 

NEWREC" . REC : =THISGRAPH . GRAPH [NEWNODE [ 2 ] ] ; 

NEWREC" . PARENT : =THISGRAPH . GRAPH [NEWNODE [ 1 ] ] ; 

NEWREC " . WT : =NEWWT ; 
.  THISGRAPH . GRAPH [NEWNODE [ 1 ] ] " . NEXTNODE 
: =NEWREC ; 
END; 
END;       {IF  NO  SELF  LOOP} 
END;       {WHILE} 
END;        {PROCEDURE  MAKEGRAPH} 


PROCEDURE  INITGRAPH (VAR  Gl:G); 
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BEGIN 

FOR  I:=l  TO  MAX  DO 
BEGIN 

NEW(G1.GRAPH[I] ) ; 
WITH  Gl . GRAPH [ I ] *  DO 
BEGIN 

NODE:=I; 
DIST:=10000000; 
PWAVE:=RANDOM; 
PDAM:=RANDOM; 
PRED : =NIL ; 
NEXTNODE:=NIL; 
DAM : =0.0 ; 
QIN:=FALSE; 
END; 
END; 
Gl. LENGTH :=0; 
END;      {PROCEDURE  INITGRAPH} 

PROCEDURE  CHOICES (VAR  Gl :G;DEST: INTEGER) ; 

TYPE  QUEUE = ARRAY [1. .MAX]  OF  BOOLEAN; 

VAR  M : INTEGER; 
Q2 : QUEUE ; 

BEGIN 

FOR  M:=l  TO  MAX  DO 

Q2 [M] :=FALSE; 
Q2 [DEST] :=TRUE; 
FOR"  M : =MAX  DOWNTO  1  DO 
BEGIN 

IF  (Q2[M]=TRUE)  THEN 
BEGIN 

IF  ( (M  MOD  10) =1)  THEN 
BEGIN 

IF  (M>9)  THEN 

Q2 [M-9] :=TRUE; 
IF  (M>10)  THEN 
Q2 [M-10] :=TRUE; 
END 
ELSE  IF((M  MOD  10)=0)  THEN 
BEGIN 

IF  (M>10)  THEN 

Q2 [M-10] :=TRUE; 
IF  (M>11)  THEN 
Q2 [M-ll] :=TRUE; 
END 
ELSE 
BEGIN 

IF  (M>9)  THEN 


Q2 [M-9] :=TRUE; 
IF  (M>10)  THEN 

Q2 [M-10] :=TRUE; 
IF  (M>11)  THEN 
Q2 [M-ll] :=TRUE; 
END; 
END; 
END; 
•FOR  M:=l  TO  MAX  DO 

Gl . GRAPH [M] A . QIN : =Q2 [M] ; 
END; 

PROCEDURE  FIND PATH ( VAR  Gl : G ; S : INTEGER ; D : INTEGER) ; 

TYPE  A=ARRAY [ 1 . . 3 ]  OF  PEDGE ; 

VAR  C: PEDGE; 

THI SNODE ,  Z , J , K : INTEGER ; 

W,TEM:REAL; 

CP:A; 

BEGIN 

Gl . GRAPH [ S ] " . PRED : =G1 . GRAPH [ S ] ; 
Gl. GRAPH [S] ~.DIST:=0.0; 
THISNODE:=S; 

FOR  Z:=l  TO  ( (D  DIV  10) -(S  DIV  10))  DO 
BEGIN 

CP [1] : =G1 .GRAPH [THI SNODE ] A .NEXTNODE; 
FOR  J: =2  TO  3  DO 

CP[J] :=CP[J-1] ".NEXTNODE ; 
" TEM:=10000000; 
K:=l; 

WHILE  ((CP[K]oNIL)  AND  (K<  =  3))  DO 
BEGIN 

W:=BUILDWTS(CP[K] ) ; 

IF  (  (CP[K]/S.REC".QIN=TRUE)  AND  (W<=TEM)  )  THEN 
BEGIN 

C:=CP[K] ; 
TEM:=W; 
END; 
K:=K+1; 
END;  . 
THI SNODE : =C" . REC" . NODE ; 
CA . REC" . DIST : =0" . PARENT" . DIST+TEM; 
C" . REC" . PRED : =C" . PARENT; 
END; 
END; 

BEGIN  {MAIN  PROGRAM} 

ASSIGN(DATAOUT, 'C:\PASCAL\OUTPUT5.TXT' ) ; 
REWRITE (DATAOUT) ; 
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RPROB ( PI , P2 ) ; 

RLATLON ( PNTLAT, PNTLON) ; 

GRAPH1. GRAPH [START] ".DIST:=0 


PDAMAGE [ 1 
PDAMAGE [ 2 
PDAMAGE [ 3 
PDAMAGE [ 4 
PDAMAGE [ 5 
PDAMAGE [ 6 
PDAMAGE [ 7 
CDAMAGE [ 1 
CDAMAGE [ 2 
CDAMAGE [ 3 
CDAMAGE [ 4 
CDAMAGE [ 5 
CDAMAGE [ 6 
CDAMAGE [ 7 
INITGRAPH 
MAKEGRAPH 


=  0 
=  0 
=  0 
=  0 
=  0 
=  0 
=  0 


858; 

21; 

087; 

0354 

0027 

0024 

0003 
=340771 
=312196 
=312196 
=129969 
=129969 
=48427 
=48427 

;graphd 
:graphd 

CHOICES ( GRAPH1 , STOP ) 
FOR  I:=l  TO  100  DO 
BEGIN 

FOR  J:=l  TO  MAX  DO 
BEGIN 

GRAPH  1.  GRAPH  [J]  ' 

GRAPHl. GRAPH [J]' 
END; 
FINDPATH (GRAPHl, START, STOP) ; 
CPNT : =GRAPH1 . GRAPH [ STOP ] " . PRED ; 
TOTAL : =FUEL (CPNT" . NODE , GRAPHl . GRAPH [ STOP] 
WRITE ( DATAOUT , STOP : 4 ) ; 

WHILE (CPNT"  .  PRED"  .NODEoCPNT"  .NODE)  DO 
BEGIN 

TOTAL : = TOTAL + 

(FUEL (CPNT" . PRED" .NODE, CPNT" .NODE) *45 
WRITE ( DATAOUT, CPNT" . NODE : 4 ) ; 
CPNT : =CPNT" . PRED ; 
END; 
WRITE ( DATAOUT , START : 4 ) ; 

WRITE (DATAOUT, GRAPHl. GRAPH [STOP] " .DIST : 10 : 2 ) ; 
WRITELN (DATAOUT, GRAPHl . GRAPH [STOP] " .DIST-TOTAL : 10 
FOR  J:=l  TO.  MAX  DO 
BEGIN 

GRAPHl.  GRAPH  [J] 
GRAPHl. GRAPH [J] 
END; 
END; 

CLOSE (DATAOUT) ; 
END. 


PRED:=NIL; 
DIST:=10000000; 


.NODE) *45.86; 


86) 


2); 


PWAVE:=RANDOM; 
PDAM : =RANDOM ; 
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INITIAL  DISTRIBUTION  LIST 


No.  Copies 
Defense  Technical  Information  Center  2 

Cameron  Station 
Alexandria  VA  22304-6145 

Library,  Code  052  2 

Naval  Postgraduate  School 
Monterey  CA  93943-5002 

Daniel  C.  Boger,  Code  OR/BO  2 

Naval  Postgraduate  School 
Monterey  CA  93943-5002 

Michael  Sovereign,  Code  OR/MS  2 

Naval  Postgraduate  School 
Monterey  CA  93943-5002 

J.  Nottenkamper  2 

ADP  Hardware  Coordinator 

Fleet  Numerical  Oceanography  Center 

Monterey  CA  93  943 

Lt  James  L.  Brown  Jr.  2 

150  Brown  Lane 
Harmony  PA  16037 

James  L .  Brown  Sr .  1 

150  Brown  Lane 
Harmony  PA  16037 
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